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We  consider  the  set  of  points  x t ^ satisfying  H(xJ  = 0, 
where  H:  -*■  is  a function  and  0 is  n regular  -'al'-'e.  This 

, t 

set,  H~‘"'(0),  is  a C"^  one-dir:ensional  nanifold,  and  each  ccnponent  can 

V 

bo  described  by  a curve  x(9).  In  this  note  a tiieorem  is  proved  •.ihich  is 

directly  related  to  and  notivatod  by  a result-^due  to  Eaves  and  Esari^ on 

piecewise  linea.-  functions.  This  theorem  relates  tie  signs  of  the  uerivatives 

X, (0)  to  the  signs  of  the  determinants  of  submatrices  of  the  Jacobian  matrix 
(1 

•H ' . Applications  to  solving  nonlinear  oouations  arc  given. 


1 . Tnt.i  oduc'' ion  and  notation 

A well  kno’vn  technique  for  solving  nonlinear  eouations  f(y)  = 0,  wiicre 
f;  -*■  R^  is  to  imbed  f into  a one-parameter  ‘’amily  of  hoirotcpy  equation: 
H(x,  t. ) = 0 where  x ^ r”  , t £ [j,  l],  H(x,  C)  H f(x),  and  n(x^,  l)  = 0 
for  sore  known  x*^  £ R^..  Examples  of  such  a function  H are 

H(x,  t}  = f(x)  - t f(x^) 

H(x,  t)  = t(x  - x^)  + (l  - t^  f(x) 


'i^e  work  of  tills  author  was  supported  in  part  by  ORR  Grant  Vo.  t 

7G-G-0'i95^in  i VCI’  Gr.ant  No.  ll.'G  76-ajOGa. 
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Many  previous  studies  have  analyzed  conditions  under  which  the  equation 
H(x,  t)  = 0 has  a solution  x(t)  which  is  a differentiable  path  for 


which  H^(x(t),  t)  is  nonsingular  for  all  t fc  [O,  l],  where 
denotes  the  derivativ*  of  H with  respect  to  the  x variables.  This 
requirement  is  equivalent  to  stating  that  the  differential  equation 

X = -H  (x,  t)“^  H (x,  t)  , x(l)  = 

X w 

has  a solution  x(t)  for  t € [O,  l],  where  x denotes  the  deri'rative 
— . Indeed,  in  this  latter  case  ^H(x(t),  t)  = 0 for  t g [O,  l]  and 
hence  H(x(l),  l)  = 0 = H(x(0),  O)  = f(x(0))  so  that  x*  = x(0)  is  a root 
of  f(x)  = 0.  One  can  then  solve  the  equations  f(x)  =0  by  integrating 
the  above  differential  equation  (assuming  x(t)  is  a differentiable  path). 
However,  the  requirement  that  x(t)  be  a differentiable  path  is  quite  severe. 

For  appropriate  discussions,  see  the  papers  of  Meyer  [9  ],  Davidenko  [1  ], 
and  Jacovlev  [ 7 ] . 

Alternative  methods  for  solving  f(x)  * 0 involve  tracking,  in  a 
limiting  sense,  points  (x,  t)  which  satisfy  H(x,  t)  = 0 by  the  so  called 

method  of  complementary  pivoting  on  a triangulation  of  R^.  All  of  these  methods 
are  extensions  of  the  seminal  work  of  Scarf  [12],  [13]  on  the  application  of 
complementary  pivoting  to  general  nonlinear  problems.  For  detailed  discussions 
see  for  example  works  of  Merrill  [8],  Garcia  [3],  and  Garcia  and  Gould  [1*],  [5l. 
as  well  as  numerous  other  papers  noted  in  these  bibliographies. 

In  general,  these  latter  complementary  pivoting  algorithms,  although 
slow  by  nature,  converge  under  assumptions  much  weaker  than  the  existence  of  a 
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diiToronf,  3 "iblf  path  x(t).  The  basic  requirenont  is  that  the  set  of  points 
(x,  t)  £ R ^ such  that  K(x,  t)  = 0 be  a onc-dir.ensional  dirf'ercn'.labi 
TT.iin  1 fold  for  which  and  x*  are  in  the  same  component,  where  f(x*)  = 0. 

Th\is  one  is  l<=“ad  to  study  differentiable  objects  such  as  x(0),  t(9),  9 g [O,  l] 

usinr;  methods  of  iiffer^.-ntial  and  conhi  natorial  tonoloF^y.  3uch  methods,  al  thouj;!’. 


hardly  new,  have  been  rarely  applied  to  the  problem  of  solving  f(x)  = 0. 

.■;oln.blo  exceptions  ai'e  the  papers  of  Saves  and  Ccarf  f2]  and  Smalc  ’■■<e  ai’e 

indebted  to  Herbert  ?carf  for  referring  us  to  the  latter  paper. 


In  this  note  a theorem  is  proved  which  lends  insight  to  the  behavior 

of  the  set  of  points  (x,  t)  for  which  H(x,  t)  = 0.  This  theorem  relates 

the  signs  of  the  derivatives  x(9),  t(9)  to  the  signs  of  the  detenriinan' s 

of  submatrices  of  the  Jacobian  matrix  H'. 

2 

Let  us  generally  consider  a C (twice  continuously  differentiable) 
function  H:  ■*  R‘^.  An  example  of  interest  is 


H(x,  t)  = f(x)  - tf(x°),  X (£  r",  t g R 


Given  y t , I*?! 


H~^(y)  = {x  € r"‘^^|k(x)  = y} 


and 


Il 


C = {x  t I rank  H'(x)  < n} 


whei-e  H'  is  the  Jacobian  rr*atrix 


I 


of  H with  respect  to  x g : 


vH  + l 


\ o t ^ 

Tile  scl  C is  said  to  be  the  set  of  criiical  nointr:  of  H,  and  H{Cj  tnc 
set  of  critical  values,  ~ H(C)  is  the  set  of  re,^u3ar  values,  oard  :> 

Tiiecrcin  fl5  ] states  that; 

Theorem  1 . Let  H:  r"  be  a map.  Then  il(C)  has 

measure  zero. 

Thus,  as  a corollary,  the  set  of  regular  values  is  dense  in  r'^  . Let  us 
henceforth  throughout  this  paper  assume  thau  0 is  a regular  value  of  H. 
The  folloving  lemma  will  be  used  [lO ] : 


,n+l 


O 

r"  be  a a map  and  let  0 be  a 


Lemma  1 . l.et  H:  R 

regular  value  of  H.  Then  H"^(0)  is  a one-dimcnsional 

manifold. 

We  now  recall  that  any  connected  one-dimensional  m.anifold  is  iiffeo- 

morphic  to  a circle  or  an  interval  (open,  closed,  or  half-open).  Thus, 

(cor.r.ected)comr-or.ent  of  ir^(O)  can  be  described  by  a curve  xCi)  which  is 

diffeomorphic  to  a cir.cle  or  an  interva]  . Furthermore,  for  any  x(0)  t H ^(0)  , 
we  have 


rank  h'  (x(0)  ) = n 


(1) 


and  x(0)  is  a unique  nonzero  vector  (x(0)  denoting  the  derivative  ol  x 


;4 
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with  respect  to  9 at  9=0).  Consequently,  we  can  differentiate 
H(x(9))  = 0 with  respect  to  9 to  obtain 

H'  (x(9))  x(9)  = 0 . ^2) 

For  a peurticular  9,  x(9)  is  a vector  tangent  to  the  curve  at  0 = 9 

and  spans  the  kernel  of  H’(x(9)). 

For  any  i = l,  2,  ...,n-t'l,  let  x^(9)  and  H^(x(9))  denote 
the  i^^  component  of  x(9)  and  the  i^*^  column  of  H'(x(9)),  respectively, 
and  let  x^(9),  H^(x(9))  be  the  remaining  components  of  x(9)  and  columns 

of  h'(x(9)),  respectively. 

2.  The  Main  Theorem 

Ouj.  following  theorem  for  maps  is  related  to  and  motivated  by  a 

theorem  of  Eaves  and  Scarf  for  piecewise  linear  maps  [2]. 

Theorem  2.  Let  H:  -►  R°  be  a map  and  0 a regular 

value  of  H.  Then  for  any  component  x(0)  of  H ^(0)  we  have 
for  all  1=1,2,  ...,n+l; 

sgn  x^(9)  * sgn  det  H^(x(9))  all  9 
or 

sgn  Xj^(0)  » -sgn  det  H^(x(0))  all  9 

(where  sgn  0^0). 


d 
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prove  olic*  theorem  in  three  parts. 


’.emirm  2.  If  x.(Q)=0  t)’.en  det  H^(x{6))  = 0. 

1 

Proof:  Py  ( 1 ) and  ( 2 ) we  have  H ' (y( 0 ) ) x( 0)  = 0 = , 

where  rank  l!'(x(6))  = n and  x(6)  0.  If  x^(0)  = 0,  then 

x^(0)  ^0  so  that  det  K^(x(0))  must  be  zero.  H 

Lerrma  3 • If  det  H^(x(0))  = 0 then  x^(9)  = 0. 

Proof:  If  det  K^(x(0))  = 0,  then  for  some  j i , we  have 
det  H'^(x(S))  ^ 0 because  by  ( 1 ) rank  H'(x(0))  = n.  For 

simplicity,  take  i = n and  j = n + 1.  Then,  x”^  ^ 

implies  x^'*’*  = H x , , . Note  that  the  last  component 

of  is  zero,  otherwise 


det  h"  = det  det  [ ] 

= det  det  [!I^,  H2,  ...,  H^]"^  [H^, 

= det  det  .■■!,  where  M 

n X n matrix.  Tlie  first 


is  the  following 
n - 1 col’jrms 


Hence,  x - 
n 


.n+1 


are  given  by  j^n-ll  = [H*^  ] 

0 j +11 

and  the  last  col'omn  is  ] ^^n+1' 

last  component  of  the  last  column  is  not  zero 

tiien  det  0,  a contradiction. 

0.  tf 
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lOTOTl'l  I*  . 


or 


x.O)  det  H^(x(0))  > 0 all  0 such  that  x.(6)  f 0 


x.(0)  det  H^(x(B))  < 0 all  0 such  that  x.(0)  ^0 


Proof:  Let  A(0)  = 


H^(x(0))  , H^(x(0))1 

^(0)^  , x.(0)  J 


L x 


B(0) 


h/(x(6))  , x^(0) 

0 , x.(0)  ^ 

where  x^(0)^  is  the  transpose  of  x^(9).  Since  x(9)  is  ortho- 
gonal to  H'(x(9)),  we  have  rank  A(0)  = n 1 for  all  0.  Since 
A(9)  is  conoinuous  in  9,  we  have  det  A(0)  >0  all  6 or 
det  A(  0 ) < 0 for  all  9 . 


Now 


det  A(9)  B(9)  = det 


[* 


L x^(9)^H^  , x(0)^  x(9)  ^ 

= x(9)^  x(0)(det  > 0 

since  x^(0)  4 0 implies  det  ^ 0 by  Lemma  3-  Thus  the 

determinants  of  A(9)  and  D(9)  have  the  same  nonzero  sign  for 
all  0 such  that  x^(0)  0.  But 

det  B(0)  = x^(9)  det  n’’(x(0)) 
which  proves  the  claim.  H 
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From  the  previous  3 lemmas,  we  get  Theorem  2 directly. 

Note  that  the  theorem  holds  if  H is  restricted  to,  say, 

R°  X [0,  l].  In  most  applications  to  x [o,  l]  a further  restriction 
would  be  required  on  the  boundaries  r’^  x {o}  and  R^  x {i} — namely,  non- 
singularity of  the  n X n submatrix  points  (x,  t)  in  the 

b Ota  da ry  for  which  H(x,  t)  * 0.  This  condition  assures  that  all  loops  which 
occur  are  contained  in  x (q^  l).  An  interesting  corollary  to  the  theorem 
is  the  following  nonotonicity  theorem. 

Corollary  Let  H:  -►  be  a map  and  0 a regular 

vadue  of  H.  Suppose  for  some  i,  ff^(x)  is  nonsingular  for  all 
X in  a particular  component  x(0)  of  H ^(0).  Then,  on  that  com- 
ponent of  H~^(0)  ,x^(9)  is  either  monotone  increasing  or  monotone 
decreasing  as  a function  of  9. 

Observe  that  ’ader  the  assumptions  of  the  Corollary  the  disting’oished  component 
of  H~^(0)  cannot  be  diffeomorphic  to  a circle. 

3.  Illustrations 
Illustration  1: 

3 2 

"Hie  theorem  is  illustrated  in  Figure  1 for  the  function  H:  R -*■  R 

given  by 

H(x,  t)  * f(x)  - tf(x°),  X 6 R^,  t g R 

where 

f^(xi.  X2)  - x^ 

fg(Xj^,  Xg)  » (x^  + Xg)^  - 10(x^  + Xg)  + 9 
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and  f(x  ) = (-1*,  9109).  T)ie  function  f^  is  a rotation  about  the  vertical 


axic  of  F(x)  = (x  - ])(x  + l)(x  - 3)(x  + 3) 


whose  graph  is 


/ _ 

I i \ v; 


' f 


w/'  r. 


BESrmiUlE  COPY 


For  this  examiiie  H (O)  is  a sint^lG  component  (x(9),  t(0)),  6 ^ [O,  l], 

—1  2 

and  the  projection  of  H (O)  into  R is  a loop.  The  matrix  H'  of 
Theorem  2 is 
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"Vf^  - f^(x°)  ' 


L**x^(xi  + x^)  - 20x^ 


- 20x, 


_Q 


9109 


According  to  the  Theorcn, 

(i)  t(9)  = 0 <“>  Vf^  and  Vf^  are  dependent  i.e.,  if  and  only  if 

3 X*  2 2 

= 0,  and  therefore  <=“>  x,  + x^  = S or  x^  = 0 . 
dx^  122 

(ii)  x^{B)  = 0 <=>  the  first  and  third  colimns  of  H'  are  dependent, 

2 2 

which  is  true  if  and  only  if  l4x^(x^  + x^)  - 20x^  = -9109/^. 

(iii)  X.  (6)  = 0 the  second  and  third  col'unns  of  H'  are  dependent, 

3^2 

which  IS  true  if  and  only  if  ^ — = 0 

dx^ 

Tne  path  shovi’n  in  Figure  1 illustrates  these  properties. 
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t - I i' ! 

Agiiin  for  Lhir  example  there  is  a sini'^lo  eom[)oncnt  in  the  projection  of  K ^(O) 

2 

into  R and  this  path  is  a loop.  It  can  be  verified  tliat  on  the  portion  of  tlio 
path  connectlnp;  to  the  "first"  sero  (point  A)  t decreases,  thon  increarer, 

then  decreases  to  0.  llcnoc  this  scpmi-nt  cannot  be  reprt'sented  as  n differentiable 
path  x(t). 


lit 


U . ApyIic"t.ions 
AppJic;ition  1 ; 

Ag  a first  arpLication  we  suTimariae  Gome  of  the  ccirjTientG  in  [ 5 ] 
on  relat  icriG  between  the  set  of  solutions  to 

H(x,  t)  = f(x)  - t f(x^)  =0,  X € t c R 

to  the  sot  of  solutions  to 


(3) 


f'(x(9)  x(9)  = - 1(9)  f(x(9)),  x(0)  = (It) 

2 0 

where  f is  a C‘~  map  with  f'(x  ) nonsinhUlar  and  > is  a real  valued 

function  of  9 satisfying  the  condition  sgn  X = sgn  det  f'(x(9)).  The 

set  of  points  satisfy-ing  equation  (3)  can  be  followed  (in  a precise  limiting 

sense)  by  a scalar  labeling  sir.plicial  pivoting  algorithm  similar  to  the 

method  introduced  by  the  authors  in  [ I*  ] . 'Hie  differential  equation  (l‘)  ’was 

introduced  by  ilmale  in  [lU]  and  has  the  credentials  of  a "global  ijewe^on  method." 

2 

Since  by  assumption  f is  C and  0 is  a reg’ular  value  of  H, 
the  solutions  to  (3)  comprise  a one-dimensisnal  manifold,  so  that  the 

component  of  ii  *(0)  containing  the  "initial  point"  (x^,  l)  may  be  described 
by  a curve 


(x(9),  t(9)),  (x(0),  t(0))  = (x  , l) 
Differentiating  (3)  wc  obtain 


f’(x(0))  x(0)  = t(0)  f(x°)  = f(x(9)) 


(5) 


V 

f: 


\ 
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if  t(G)  ^ 0.  ThuG  if  t(0)  ^0  it  followr.  from  Theorem  2 that 


det  f(x(0))  / 0 so  that 


x(e)  = Ijll  f'{x(e))-^  f(x(0))  . 


In  this  particular  case  cur  theorem  says 


sgn  t(9)  = sgn  det  f'(x(0))  all  9 


sgn  t(0)  = -sgn  det  f'(x(9))  all  9 


We  now  observe  that  (5)  provides  a special  instance  of  (h)  if 
t ( 0 ) 

17'o \ = -sgn  det  f'(x(6)).  Recall  that  t(0)  = 1,  and  adopt  the  ccnventicn  that 
1 1 0 ; 

if  det  f'(x^)  > 0 (<  O)  we  initially  move  on  the  path  in  such  a iirectior.  that  t 

“t  ( 0 ) 0 

decreases  (increases)  from  its  initial  value  1.  Then  sgn  -7^-r  = -sgn  det  f'(x  ) 

u V U J 

and  by  Theorem  2 this  will  be  true  for  all  0 such  that  t(G)  > 0.  This 
proves  that  a piece  of  the  solution  to  (3)  is  a solution  to  (I4),  for  Eralc's 
algorithm  tercninates  the  moment  the  first  zero  is  encountered  (t(G)  = 0) 
whereas  (3)  c;an  be  "continued."  In  Figure  1,  the  path  from  x^  to  the  point 
labeled  A is  the  Sraaie  path  (solution  to  (U)),  The  solution  to  (3)  is  the 
entire  closed  loop  which  contains  all  of  the  zeros  of  f. 


Now  consider  Figures  3 and  . Here  we  are  solving  the  equation: 


f^(x^,  x^)  = - 3 = 0 

f^Cx^,  x^)  = (x^  f - 10(x^  + x^)  - 3 = 0 . 


l6 


The  obvious  relation  to  the  function  f in  Illustration  1 should  be 
noted.  It  is  neon  in  Figure  3 that  the  projection  of  H ^(O)  into 
R contain.s  two  components.  That  component  pa.ssin/’;  throur.h  the  "starting 
point"  docs  not  pass  throuph  a root.  However,  consider  Figure  U where 

the  starting  point  is  chosen  "at  infinity."  In  tJiis  case  the  projection 

of  H (O)  into  R is  a single  path  x(0)  which  passes  through  both 

roots  of  f.  It  can  be  verified  that  this  interesting  behavior  will  occur, 
for  this  example,  with  any  x^  sufficiently  large  in  norm. 

In  [ll<]  Smale  has  demonstrated  a result  which  can  be  restated  as  follows 

Supuose  there  is  a bounded  open  set  C such  that  f;  C -*•  r", 


c 

and  3C  are 

connected. 

3C 

smooth , and  x G 

3C  =»>  det  f ' ( X ) 

> 0 

(f 

’(x))"^  f(x) 

intersects 

3C 

transvcrsally  at 

X.  Suppose  f 

€ C' 

x*^  £ 3r  and  0 is  a regular  value  of  H(x,  t)  = f(x)  - t f(x^).  Then 
the  connected  component  of  H ^(0)  containing  (x^,  l)  will  contain  a zero 
of  f. 

Figure  U demonstrates  a set  C satisfying  Gmale's  assumptions.  Note 


that  only  1 of  the  roots  is  contained  in  this  set. 
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Application  2; 

In  [6  1,  Ci.arcia  and  Zan,<wilJ  r.howed  how  to  find  all  solutions  to 
certain  systems  of  n nonlinear  equations  in  n complex  variables.  For 
example,  they  show  how  to  construct  all  solutions  to  an  arbitrary  system  of 
n polynomial  equations  in  n unknowns. 

The  underlying  theorem  in  [ 6 ] proves  a mcnotoriic  behavior  of  the  paths 
of  solutions  to  a particular  set  of  homotopy  equations.  This  result  is 
achieved  by  use  of  the  Cauchy-Reimann  conditions  to  show  that  det  H^(x,  t)  > 0 
for  any  (x,  t)  in  ^ [O,  l],  (where  is  the  Jacobian  of  H written 

as  a function  in  x [O,  l]  The  result  is  therefore  a special  case  of 

our  corollary  to  ITieorem  2,  which  in  this  instance  states  that  the  path  must 
be  monotonic  nondecreasing  (or  monotonic  nonincreasing)  in  the  variable  t. 

This  key  result  may  then  be  used  in  [6]  to  show  that  starting  from  any  solution 
to  H(x,  0)  = 0,  the  path  x(t)  satisfying 

H(  x(t ) , t ) = 0 

must  yield  an  x(l)  which  solves  the  given  system  of  equations.  See  [6  ] for 
a complete  treatment  of  this  problem. 

Application  3: 

Let  us  consider  a C function  f;  r"^  -*■  where  det  f'(x)  >0 
all  x and  lim  ||f(x)|l  = “>  (f  satisfying  the  latter  condition  is  said 

ilxll- 

to  be  norm-cperci ye ) . 

This  condition  on  f is  essentially  that  for  the  lladamard  theorem  (see 
Theorems  'j.'i.O  and  9. 3*10  of  [H].  See  also  Theorem  10. U. 3).  It  is  known  tluit 
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for  .my  y £ , the  above  aosumptiono  ar.oure  the  existence  of  a unique 

X*  sat 

f(x*)  = y . 

To  find  X* , consider  the  homotopy 

H(x,  t)  = f(x)  - [t  f(x°)  + (1  - t)  y]  = 0 (x  € R",  t 6 R) 

lor  ail  arbitrary  x^  6 R^  (This  is  the  limiting  path  generated  by  the 
cor.plenentar:»'  pivoting  algoritlim  of  Garcia-Gould  in  [ ^]).  TYien 

det  H^(x,  t)  = det  f'(x)  > 0 

for  any  (x,  t),  so  that  the  algorithuj  of  Garcia-Gould  will  trace  a ' 

path  (x(9),  t(9))  which  by  the  corollary  to  Theorem  2 will  be  nonctonic 

in  t.  This  monotonicity,  along  with  the  assumption  lim  |(f(x)||  = 

I |xl  l-«° 

implies  that  t cannot  be  asymptotic.  One  is  therefore  assured  of  finding 
the  unique  solution  (x(0),  t(0))  = (x*,  O)  of 

f(x)  = y 

for  an  arbitrary  y € r'*. 
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